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ABSTRACT 

We compare DNS calculations of homogeneous isot ropic turbulen ce with the statistical properties 
of intra-cluster turbulence from the Matryoshka Run (IMiniatill2014fl and find remarkable similarities 
between their inertial ranges. This allowed us to use the time dependent statistical properties of intra¬ 
cluster turbulence to evaluate dynamo action in the intra-cluster medium, based on earlier r esults from 
numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo (|Beresnvakl2012ll . We argue 
that this approach is necessary (a) to properly normalize dynamo action to the available intra-cluster 
turbulent energy and (b) to overcome the limitations of low Re affecting current numerical models 
of the intra-cluster medium. We find that while the properties of intra-cluster magnetic field are 
largely insensitive to the value and origin of the seed field, the resulting values for the Alfven speed 
and the outer scale of the magnetic field are consistent with current observational estimates, basically 
confirming the idea that magnetic field in today’s galaxy clusters is a record of its past turbulent 
activity. 

Subject headings: cosmology: theory—magnetohydrodynamics—MHD dynamo 


1. INTRODUCTION 

The hot intracluster medium (ICM) of galaxy clusters 
(GO) is well known to be magnetized from radio obser¬ 
vations. These reveal both the occurrence of Faraday 
rotation effect on polarized radiation fro m background 
quasars (jClarke et al.l j2001|: |Glarkf jl200^ and of diffuse 
synchrotron emission ( Ferrari et al.ll200^ from the ICM. 
Estimates of the magnetic field based on these observa¬ 
tions range between a fraction and several fiG. Measure¬ 
ments on the structural and spectral features are sparse 
and more difficult, but indicate steep power-laws be¬ 
low f ew tens of kpc (|Laing et al.ll200^ iKuchar fc EnfilinI 
[Mill . For massive clusters, turbule nce in the ICM is 
main l y driven by structure formation (iN orman fc Brva3 
1999 : iRvu et ail [20081: iVazza et al.l 120111 : iMiniatil 12014 
201511 . The most important magnetic field amplifica¬ 


tion mechanism in the ICM is the small scale or fluc¬ 
tuation dynamo (SSD), operating on scales smaller than 
the turbulence outer scale. Kinematic regime of SSD, 
i.e. when the back reaction of the magnetic field on the 
flow is negligible, has been studied in great detail pre¬ 
viously (lKazantsevlll968t iKraichnan fc NagaraianlllOOTt 
iKulsrud fc AndersonI Il992| l. In kinematic regime the 
magnetic energy grows exponentially, till the approxi¬ 
mation breaks down, roughly in a dynamical time multi¬ 
plied by where Re is an effective Reynolds num¬ 

ber. The extremely hot and rarefied plasma of the cluster 
have very large collisional mean free paths, around 


A ~ 10^ pc(n/3 X 10 


-3cm-3^-l 


)-i(r/10keV)3/2^ (1) 


at the same time, given the observable magnetic fields 
around 3 /xG, the Larmor radius is smaller by many or¬ 
ders of magnitudes: 


CL « 10-®pc(r/10keV)(R/3/xG)-L (2) 


Such situation, known as “collisionless plasma” is 
challenging from theoretical viewpoint, since nonlinear 
plasma effects are dominating the transport, which has 
been known since early Lab plasma experiments, when 
it became clear that collisional “classic transport” is 
gross ly insufficient to explain cross field diffusion (see, 
e.g., iCaleev fc SagdeevI Il979fl . As a rule of thumb, 
the actual effectiye parallel mean free path is smaller 
than the one obtained by collisional formula, but larger 
than the Bohm estimate (Ae// ~ r^). The search for 
this “mesoscale” for cluster conditions resulted in esti¬ 
mates for the mean free path of the proton in th e ICM 
around 10~^ — I0~^pc (fScheko chihin fc Cowle^J 2006; 
Beresnyak fc Lazarianl | 2006l : LSchekochihin et al.f 2008; 
Brunetti fc LazariaDll20lfh . From these estimates we ex¬ 
pect clusters to be turbulent with Reynolds numbers Re 
exceeding 10^^. Combining this with the aboye estimate 
of the kinematic SSD growth rat es, for a dynam ical time 
~ eddy turnoyer time ~ 1 Gyr (|Miniatil 120141 1. we esti¬ 
mate that the exponentiation timescale will be smaller 
than 1 Gyr « 1 kyr. 

The remainder of this paper is organized as follows: in 
Section 2 we discuss the properties of nonlinear regime 
of the small-scale dynamo which is supposed to domi¬ 
nate during most of the cluster lifetime; in Section 3 we 
point to the inadequacy of current MHD cosmological 
simulations, as far as dynamo is concerned, and suggest 
a different approach; in Section 4 we describe new homo¬ 
geneous dynamo simulations with intermittent driying; 
in Section 5 we explain our cosmological hydrodynamic 
model of the cluster; in Section 6 we combine the knowl¬ 
edge obtained in previous sections and analyze cluster 
simulations to derive the properties of the cluster mag¬ 
netic fields; in Section 7 we discuss implications and com¬ 
pare with previous work. 
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2. NONLINEAR SMALL-SCALE DYNAMO 

As the kinematic approximation of SSD breaks down 
very quickly, the dynamo spends most of the time 
in the nonlinear regime. In this regime, inclusive of 
the back reaction of the magnetic field on the flow, 
the magnetic energy continues to grow as it reaches 
equipartition with the turbulent kinetic energy cas¬ 
cade at progressively larger scales (jSchliiter fc BiermannI 
Il95flfl . At this stage the magnetic energy is char¬ 
acterized by a steep spectrum and an outer scale, 
LB , a small fraction of the kinetic energy outer 
scale (iHaueen et al.l 12004 iBrand enb uxg fc Subramanianl 
l200,^l^n et al.ll2008l : lCho et aLll200^ . This picture has 
been later argued to be true in any high-Re flow, with 
the ar gument relying o n locality of energy transfer func¬ 
tions (jBeresnvakll20I2t l. It also followed from this study 
that the growth rate of the magnetic energy corresponds 
to a certain fraction of the turbulent dissipation rate, 
with this fraction being a universal dimensionless num¬ 
ber around 0.05, and tat the m agnetic outer sc ale Lb 
grows with time as Lb ^ (|Beresnvakl l2012h . The 
growth of magnetic energy reaches final saturation when 
Lb IS a substantial fraction of the outer scale of the tur¬ 
bulence. However, this never happens in clusters, as we 
show below. 

3. LIMITATIONS OF COSMOLOGICAL DYNAMO 
SIMULATIONS 

An important implication of the above picture is that 
the memory of the initial seed field is quickly lost and the 
cluster magnetic field is expected to depend only on the 
cluster turbulent history. While this theoretical insight 
was certainly useful, its applications to cluster formation 
were not immediately realized. There are two main rea¬ 
sons for this. Firstly, while there has been considerable 
progress in computational models of structure formation, 
and GCs in particular, the level of dynamic range of 
spacial scales achieved so far is considerably below the 
threshold necessary for the turbulent dynamo to oper¬ 
ate efficiently. In fact, numerical MHD models of GCs 
typically report rather weak magnetic field amplification 
roughly by factors < 30 j Miniati et al.ll2001l :lDolag et al.l 

l200^[Dirbois fc Tevssieii2008HVazza et al.ll2014l l. includ- 

ing significant contribution from adiabatic compression. 
As alluded aboye, the reason is ascribed to the low Re 
of the simulated flows. The kinematic growth rate is 
q « Re^^^/SOrr (iHai^en et al.ll2004lSchekochihin et al.l 
12004 lBeresnyakll2012f l. where rr is the turnoyer time of 
the largest eddy. So eyen with Re ~ seyeral xIO^, typi¬ 
cal for cluster simulations, the dynamo will be stuck for 
seyeral dynamical times in a kinematic regime, i.e. sey¬ 
eral Gyr, while in nature this stage will be many orders of 
magnitude quicker than the dynamical time (see also Sec¬ 
tion I). Fig. [T] demonstrates the difference between the 
magnetic energy growth between the case with yery large 
Re (straight line) and Re that are ayailable with current 
numerical capabilities (actual growth obtained in simu¬ 
lations with Re=1000 and 3300). The growth observed 
in simulations is delayed due to the grossly prolonged 
kinematic stage. Secondly, in view of the current under¬ 
standing of MHD dynamo (Section [5]) , lack of detailed 
knowledge about the IGM turbulence precludes accurate 
estimates of both the magnetic energy and, in particular. 



Fig. 1.— The growth of dynamo can be significantly delayed 
in simulations compared with actual high-Re flows, which is de¬ 
picted here. Here we compare the growth of magnetic energy from 
4 X 10“^ below equi partition in Re = 10® and Re = 6.6 X 10® simu¬ 
lations (solid) from iBeresnyaM II2012I1 and the hypothetical Re=oo 
case which is represented by linear growth Cset (dashed). If we 
stop simulation around the dashed line, which roughly correspond 
to the duration of cosmological simulation ~ 11 dynamical times, 
the magnetic energy will be grossly underestimated and its value 
will depend on both the seed field and the effective Reynolds num¬ 
ber of the simulation. 

the outer scale of the magnetic spectrum. 

Below we report on the progress with the approach 
which is different from direct approach of cosmological 
MHD simulation, which, given present state of our nu¬ 
merical capabilities, as we argued above is completely 
inadequate. We have recently employed a novel tech¬ 
nique to model the formation of a massive GG wi th suffi¬ 
cient resolu tion to resolve the turbulent cascade (iMiniatil 
12014 I2015H . We have extracted the time dependent 
properties of the turbulence and used this information in 
combination with independent results on turbulent dy¬ 
namo obtained from high resolution periodic box simu¬ 
lations. The novelty and advantage of our approach is 
that the turbulence is self-consistently estimated through 
a numerical hydrodynamic model of structure formation, 
while the magnetic field evolution is estimated based on 
theory, which was confirmed in large-scale homogeneous 
dynamo simulations, robustly tested by studying low Re 
effects in a scaling study. Importantly enough, such dy¬ 
namo simulations, unlike cosmological cluster models, 
are not limited in the number of dynamical times one 
can simulate. 

4. DYNAMO SIMULATIONS 

We have extended the study of statistically ho- 
m ogeneous i s otrop ic small-scale dynamo simulations 
in iBeresnv^ (I2012D with a series of simulations with in¬ 
termittent energy injection into the velocity field, with 
the period 1,2,4, and 8 self-correlation timescales of ve¬ 
locity, Tc- All simulations have magnetic Prandtl num¬ 
ber Prm = 1 and driving in Fourier space was limited 
to lower harmonics (|fc| < 2.5). We started each MHD 
simulation by seeding low level white noise magnetic field 
into the dataset obtained from driven hydrodynamic sim¬ 
ulation which reached statistically stationary state. This 
dataset was further evolved by full incompressible MHD 
equations. Figure [2] shows the evolution of magnetic 
energy in time. The previously measured normalized 
growth rate Ce = 0.05 is roughly consi stent with most o f 
the data. An important prediction of iBeresnvakI (120121) 
was also that the magnetic outer scale is proportional to 
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Fig. 2. — We plotted the magnetic energy, jl (Alfven units), 
divided by the energy driving rate, e, resulting in I2e (simu¬ 
lation time units) versus time (sim ulation time units), similar to 
results reported in IBeresnvald 1)20121') . Several curves correspond to 
intermittent driving with the period of driving varied from 1 to 8 
Tc and the half of the period the driving was on while the other 
half off. We also put a simulation with constant driving for com¬ 
parison. Two thin lines correspond to the conversion efficiencies of 
nonlinear dynamo Ce of 0.045 and 0.05 (dimensionless). 

v\/e and grows in time as Since we are going to 

use this conjecture to estimate the outer scale of cluster 
magnetic fields, we plotted the v\/e versus the magnetic 
outer scale, which we determined from the peak of mag¬ 
netic spectrum. The constant driving simulation (upper 
panel of Figure [3]) showed good agreement with the pro¬ 
posed scaling and we have determined the dimensionless 
coefficient ci in the relation Lb = civ^je to be around 
0.2 (best fit 0.18). The intermittently driven simulation 
have shown large scatter which is due to the fact that 
turbulence spectra do not depend instantaneously on the 
energy injection rate, but have a memory of the previ¬ 
ous state around about one dynamical time. Also, the 
cascade rate at the equipartition scale is delayed com¬ 
pared to the injection rate. We found that averaging 
cascade rate over 2tc and introducing phase delay of 7 r /2 
will work the best to reproduce the Lb relation and we 
presented the plot of such constructed Lb - v\/e on the 
lower panel of Figure [3l The scatter was significantly 
reduced and the derived c; coefficient is also the same as 
for the constant driving case. 

5. CLUSTER SIMULATIONS 

We use the Matryoshka run to extract the time de¬ 
pendent turbulence properties of the ICM of a mas¬ 
sive GC, with total mass at redshift z = 0 of 

1.3 X lO^^Mrr^, forming in a concordance A-CDM uni¬ 
verse (iKomatsu et al.ir2009f) . The simulation was carried 
out with CH ARM, an Adaptive-Mesh- Refinement cosmo¬ 
logical code (iMiniati fc Colella 1200711 . We use a concor¬ 
dance A-CDM universe with normalized (in units of the 
critical value) total mass density, Llm = 0.2792, bary- 
onic mass density, Of, = 0.0462, vacuum energy density, 
Da = 1 — Lljn = 0.7208, normalized Hubble constant 
h = iJo/IOO km s“^ Mpc“^ = 0.701, spectral index of 
primordial perturbation, = 0.96, and rms linear den¬ 
sity fluctuation within a sp here with a comoving radius 
of 8 h~^ Mpc, (Ts = 0.817 (jKomatsu et al.l [2009( 1. The 
simulated volume has comoving size of Lbox = 240 
Mpc on a side. The initial conditions were generated 



Fig. 3. — The relation between magnetic energy and the outer 
scale in the DNS of nonlinear small scale dynamo. The upper plot 
correspond to th e case with constant driving, reported earlier in 
IBeresnvakI II2012I') . while the lower plot has been produced in oth¬ 
erwise similar simulation with intermittent driving with period of 
Stc (self-correlation timescales). In the case of intermittent driving 
we used the dissipation rate e averaged over 2tc in a manner similar 
to the analysis of cluster data. We defined outer scale of magnetic 
field through its peak wavenumber of the spectrum Lb = 27r/fcmax • 
The best fit corresponds to the coefficient ci 0.18 in Eq. {5jl- 

on three refinement levels with graf ic++ (made publicly 
available by D. Potter). For the coarsest level we use 
512^ comoving cells, corresponding to a nominal spatial 
resolution of 468.75 h~^ comoving kpc and 512^ particles 
of mass 6.7 x 10® Mq to represent the collisionless 
dark matter component. The additional levels allow for 
refined initial conditions in the volume where the galaxy 
cluster forms. The refinement ratio for both levels is, 
nfgf = Axi/Ax £+1 = 2, £ = 0,1. Each refined level cov¬ 
ers 1/8 of the volume of the next coarser level with a 
uniform grid of 512^ comoving cells while the dark mat¬ 
ter is represented with 512^ particles. At the finest level 
the spatial resolution is Ax = 117.2 h~^ comoving kpc 
and the particle mass is 10® h~^ M 0 . As the Lagrangian 
volume of the galaxy cluster shrinks under self-gravity, 
three additional uniform grids covering 1 /8 of the volume 
of the next coarser level were employed with 512®, 1024® 
and 1024® comoving cells, respectively, and nfgf = 2,4, 2, 
for £ = 2,3,4, respectively. All of them were in place by 
redshift 1.4, providing a spatial resolution of 7.3 h“^ co¬ 
moving kpc in a region of 7.5 h“® Mpc, accommodating 
the whole virial volume of the GC. The ensuing dynamic 
range of resolved spatial scales is sufficiently large for the 
emergence of turbulence. The resul ts of the cluster sim - 
ulation is described in full detail in (IMiniatill20l412015fl . 

6 . ANALYSIS OE CLUSTER SIMULATIONS 

Using a Hodge-Helmholtz decomposition it was found 
that between 60 and 90% of the kinetic energy of the clus¬ 
ter t urbulence is in the soleno i dal c omponent (|Miniatil 
120151 see also iFederrath et all (|2011IH . This is the rel¬ 
evant component for the discussed small-scale dynamo 
mechanism and the key question is whether it resembles 
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distance, Mpc 


Fig. 4. — We used structure functions to test statistical isotropy 
of the given cluster simulations (upper panel). Also we compare the 
relation between second order total structure function and the third 
order parallel structure functions in the cluster simulation (solid 
line) and in a periodic box simulation of statistically stationary 
driven turbulence (dashed line), after rescaling the dissipation scale 
to the same number (lower panel). Note that no fitting has been 
involved in making the both figures. We conclude that the second 
order structure function in the range 0.14-0.4 Mpc could be used 
to robustly estimate the turbulent dissipation rate associated with 
the solenoidal velocity component. 


homogeneous isotropic turbulence in the inertial range. 

In Fig m upper panel we checked statistical isotropy 
of the cluster turbulence. We compared the longitudi¬ 
nal velocity structure function (SF) with the analytical 
expression that presumes statistical isotropy. Statisti¬ 
cal isotropy seems to be satisfied quit e well on all sc ales 
of interest consistent with results in (|Miniatill2015D . A 
more critical test is provided by the relation between 
structure functions of different order. For example the 
dimensionless ratio / {{Sv\\i)^) is of interest to 

relate the energy cascade rate with the energy content 
of the cascade. In the lower panel of Fig 0] we stud¬ 
ied the comparison between this ratio in the cluster sim¬ 
ulation and in the homogeneous incompressible driven 
turbulence. For the latter we used data from fully re¬ 
solved direct numerical simulation of incompressible hy- 
dro dynamic driven turbul e nce i n a periodic box, see, 
e.g. iBeresnvak fc LazarianI (j2009| l. Both cluster and box 
simulation exhibited a clear well-pronounced dissipation 
interval which we used to convert box simulation units 
into physical scale units of the cluster simulation. Note, 
however, that no fitting has been involved on the y-axis. 
Given relatively short inertial range the correspondence 
between homogeneous isotropic turbulence statistics and 
cluster statistics is quite remarkable. From the above 
comparison we conclude that the second order structure 
function of the cluster simulations in the range of scales 
0.14-0.4 Mpc could be reliably used to estimate the tur¬ 
bulent dissipation rate, eturb, associated with the incom¬ 
pressible velocity component and necessary to evaluate 
dynamo action in the ICM. 

The turbulence dissipation rate is then estimated as 
follows: 

eturb = (ci/c2)(5/4)((5ui)^)^/^/Z, (3) 

where C 2 ~ 27 is the ratio of the structure functions 
reported in Fig. 0 ] and ci ~ 1.17 is a factor to correct 
for dissipation effects, as in our finite Re simulations the 


redshift 



time [Gyr] 

Fig. 5. — The evolution of the turbulent dissipation rate (upper 
plot), Alfven speed (middle) and the outer scale of the magnetic 
field (bottom) inferred from the cosmological cluster simulation by 
using self-similar laws for turbulence and dynamo tested and nor¬ 
malized in well-resolved high-ite DNS. In obtaining we divided 
the magnetic energy Eg from Eq. [J] to the average density within 
1/3 of the virial radius. 

Kolmogorov’s -4/5 normalization slightly underestimates 
the turbulent dissipation rate. 

As expected, at a given time tturb is a rather constant 
function of I within the inertial range. The observed 
deviation was used to estimate the error of the measure¬ 
ment of eturb- We plotted the dissipation rate determined 
in this manner on the top panel of Fig [5l We used the 
velocity structure function calculated within 1/3 of the 
virial radius of the simulated cluster for each data-cube. 
As we see from this figure, the dissipation rate varies 
non monotonically over roughly an order of magnitude 
in scale over the lifetime of the cluster. The errorbars 
defined above indicate the deviation from Kolmogorov’s 
self-similarity and were rather small, except for the time 
intervals where the rate was changing rapidly, i.e. the 
cluster was either relaxing of experiencing a fresh injec¬ 
tion of kinetic energy. 

We then estim ated magnetic energy density 

as (jBeresnvakll2012D 

— / CEpeturbeli (4) 

Jo 

with Ce = 0.05. We plotted Alfven velocity on the 
middle panel of Fig [5l Furthermore, as was shown in 
iBeresnvM (|2012f ). for statistically stationary turbulence 
the magnetic energy containing scale could be estimated 
as 

Lb — etiVj^j Cturb^ (b) 

where ci ~ 0.18 is a universal coefficient, which could be 
determined in DNS, see Fig. 0 ] Our cluster turbulence 
was rather non-stationary, however, as discussed in Sec¬ 
tion 4, the estimate Eq. (01) can also be applied to non- 
stationary driven turbulence as long as the dissipation 
rate is averaged over a timescale around one dynamical 
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time, see below. This is because hydrodynamic cascade 
has a memory over around one dynamical time and the 
changes in the driving rate do not instantaneously affect 
turbulent rate on small scales (Section 4). So, in using 
Eq. ([U we used the eturb averaged over 2 Gyr, which 
approximately corresponds to two dynamical times. 

The middle and bottom panels of Fig [S] show time 
evolution of the average RMS Alfven speed, va = 
{2EbI and the magnetic outer scale Lb- Note that 
while Va grows monotonically, Lb can decrease some¬ 
what during prolonged increase of the turbulent activity, 
such as during several major mergers. 

Our estimates for z ^ 0 characteristic values of 
Alfvenic speed va ~ lO'^cm/s and the outer scales 
Lb '^ 30 — 50 kpc, are consis tent with the observed values 
reported in the literature (iMcNamara fc NulsenI 120071: 
Eilek fc Owenl[2(io3 iBonafede et al.l 120101: iGovoni et al l 
2 nn(Ti 1201 mi ." This i ndicates that the type of nonlinear 
dynamo described in iBeresnveiS (j2012l l is probably oper¬ 
ating in clusters, while kinematic models would be chal¬ 
lenged to achieve this. 

One interesting conclusion from our results on Fig [5] 
is that the outer scale of the magnetic field grows rel¬ 
atively quickly after the beginning of the simulation. 
This is different from direct MHD cluster simulations 
that have mostly kinematic growth with a magnetic 
spectrum peake d on numerical dissipation scale, e.g., 
iXu et all (1201211 . Note that the scale of the magnetic 
field plays crucial role in cosmic ray escape times, there¬ 
fore correctly estimating magnetic outer scale is essen- 
tial for models of particle a cceleration in cluster s (see , 
Brunetti fc LazarianI l200'^ iBrunetti fc LazarianI 120111 : 
Beresnvak et al.ll2013t lMiniatill2015fl . 

7. DISCUSSION 

Similar idea based on post -processing of hyd rodynamic 
data was also employed in iRvu et aTT (j2008ll , but with 
substantial differences. The turbulence in these early cal¬ 
culations was not as resolved as in ours, and the growth 
of magnetic energy and Alfven scale were not estimated 
from the turbulent dissipation rate and the precise esti¬ 
mate of Cb, as we did here. 

One of the differences between cluster turbulence and 
the kind o f statis tically stationary turbulence studied in 
iBeresnv^ (j20I2t l was the strong variations of the cas¬ 
cade rate over timescales of 1-2 dynamical timescales of 
the cluster. Our estimates of the efficiency in the case 
of intermittent driving from this work are roughly com¬ 
patible with Ce = 0.05 and further work with higher 
Re is expected to clarify whether the differences between 
constant and intermitted driving are significant. We con¬ 
cluded that the effects of intermittent driving could prob¬ 
ably be ignored at the level of precision of the eturb mea¬ 
surement. 

The actual calculation for the evolution of va and Lb 
was started at time 4.5 Gyr. This artificially assumes 
that eturb was zero for all times earlier than 4.5 Gyr. 


However, we find that despite this fairly unrealistic as¬ 
sumption, the values of va and Lb quickly converge to 
the asymptotic values and, as we argued above, this ini¬ 
tial state is quickly forgotten. All basic properties of 
the cluster, such as its mass, size and thermal energy 
continue to grow along with its magnetic energy and 
magnetic outer scale. The detailed comparison between 
thermal, turbulent and magnetic energy components of 
the cluster has been performed in our companion pa¬ 
per (|Miniati fc Beresnvakll20i^ . There it is found that 
the fraction of the thermal energy arising from the tur¬ 
bulent dissipation rate changes relatively little over the 
cosmological time and the turbulent Mach number is also 
rather stable. Since the magnetic energy is also a fraction 
of the accumulated turbulent dissipation rate, the plasma 
P in our cluster fluct uates around a constant va lue ^ 40 
for the past 10 Gyr (jMiniati fc Beresnvakll2015ll . 

Our treatment of cluster turbulence with ILES, as well 
modeling th e evolution of the magnetic energy with the 
model from iBeresnv^ (|2012f) relies on an assumption 
that the Reynolds numbers in clusters are high. For ex¬ 
ample, our comparison of the cluster simulation and the 
DNS leads to an estimate of an effective Kolmogorov (dis¬ 
sipation) scale for the cluster simulation of rj Ki 2.7 kpc, 
corresponding to an effective Re around 3000. We actu¬ 
ally expect clusters to have higher Re, as briefly discussed 
in Section 1, due to the collect ive microscopic scatter¬ 
ing in the high-<3 ICM plasm a (iSchBbjchjMiLJgAlQwley 
2006; iLazarian fc Beresnvakl [20061 iSchekochihin et al. 


200S ; IBrunetti fc LazarianI l20 1 ifT ^An important obser¬ 


vational test to the problem of the IGM viscosity are the 
measurements of Faraday rotation in AGN sources lo¬ 
cated in clusters, which allowed to probe sub-kiloparsec 
scales due to relatively high resolution of radio maps 
(iKuchar fc Enfilinll2011l : iLaing et'nil2008l: IGovoni et al.l 
12010(1 . The inferred magnetic spectrum in these mea¬ 
surements is negative and steep, typically around Kol¬ 
mogorov in the range of scales below 5 kpc and down 
to the resolution limit. Such a magnetic spectrum is 
expected from MHD turbulence with small dissipation 
scales. It would be grossly inconsistent with magnetic 
spectra obtained in either kinematic dynamo models, due 
to their positive spectral indexes, or with MHD mod¬ 
els using Spitzer viscosity, which would typically give 
rather sha ll ow sp ectrum with index around —1, see, e.g. 
iGho et ahl (|2002tl . We conclude that even though it is 
quite obvious that magnetic Prandtl numbers in the IGM 
are very high, the viscosity is not large enough to affect 
magnetic spectrum above 1 kpc. Therefore, just from 
this observational constraint we expect the Re in clus¬ 
ters to be at least 10^ and probably much higher. Our 
calculation relied on this fact and the results, grossly 
consistent with the current observational properties of 
clusters, provide another support for the picture of a tur¬ 
bulent ICM, as opposed to an earlier view of a viscous 
and laminar IGM. 
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